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ABSTRACT 

We have carried out a deep imaging survey for Ha emitting galaxies at z ~0.24 using a narrowband filter tuned 
with the redshifted line. The total sky area covered is 706 arcmin 2 within a redshift range from 0.234 to 0.252 
(&=0.018). This corresponds to a volume of 3.9 x 10 3 Mpc 3 when O matter = 0.3, fl\ = 0.7, and Hq=70 km s -1 
Mpc" 1 are adopted. We obtain a sample of 348 Ha emitting galaxies whose observed emission-line equivalent 
widths are greater than 12 A. We find an extinction-corrected Ha luminosity density of 10 39,65 -° 12 ergs s" 1 Mpc" 3 . 
Using the Kennicutt relation between the Ha luminosity and star formation rate, the star formation rate density in 

+<'.(>!)h , , , .,.-1 Mp C -3 This value is higher by a factor of 3 than the local 



the covered volume is estimated as 0.036^{j ^ Mq yr 
SFR density. 

Subject headings: galaxies: distances and redshifts — 
function 

1. INTRODUCTION 

The star formation rate density is one of the important ob- 
servables for our understanding of galaxy formation and evolu- 
tion. Madau et al. (1996) compiled the evolution of star for- 
mation rate (SFR) density, psfr, as a function of redshift from 
the local universe to high-redshift. It is known that psfr steeply 
increases from z — to z ~ 1 (Gallego et al. 1995; Tresse & 
Maddox 1998, hereafter TM98; Songaila et al. 1994; Ellis et 
al. 1996; Lilly et al. 1996; Hogg et al. 1998). However, p SPR 
appears to be flat beyond z = 2 (Madau, Pozzetti, & Dickin- 
son 1998; Pettini et al 1998; Steidel et al. 1999; Barger et al. 
2000; Fujita et al. 2003a; see also Trentham, Blain, & Goldader 
1999). 

Although recent investigations have been devoted to new es- 
timates of psfr at high redshift, it seems also important to pro- 
vide new estimates in nearby universe because only several in- 
vestigations have been available between z — to z ~ 1 . For 
example, psfr at z — 0.2-0.3 is based solely on the observa- 
tion by TM98. In this Letter, we present a new estimate of psfr 
at z ~ 0.24 based on the largest sample (more than 300 galax- 
ies) studied so far, which was obtained with Suprime-Cam on 
the Subaru Telescope. 

Throughout this paper, magnitudes are given in the AB sys- 
tem. We adopt a flat universe with fi ma tter = 0.3, £7a = 0.7, and 
h = 0.7 where h = H /(100 km s" 1 Mpc" 1 ). 

2. OBSERVATIONS AND DATA REDUCTION 



We have made a very deep optical imaging survey of Ha 
emitters at z ~ 0.24 in a sky area surrounding the SDSSp 
J104433.04-012502.2 at redshift of 5.74 (Fan et al. 2000; Djor- 
govski et al. 2001; Goodrich et al. 2001) using the Suprime- 
Cam (Miyazaki et al. 2002) on the 8.2 m Subaru telescope 

1 Based on data collected at Subaru Telescope, which is operated by the National Astronomical Observatory of Japan 

2 Astronomical Institute, Graduate School of Science, Tohoku University, Aramaki, Aoba, Sendai 980-8578, Japan 

3 National Astronomical Observatory, Mitaka, Tokyo 181-8588, Japan 

4 Department of Astronomy, Graduate School of Science, University of Tokyo, Tokyo 113-0033, Japan 

5 Research Center for the Early Universe, School of Science, University of Tokyo, Tokyo 1 13-0033, Japan 

6 Subaru Telescope, National Astronomical Observatory, 650 N.A'ohoku Place, Hilo, HI 96720 

7 This survey is originally intended to search for Lyo emitters at z ~ 5.7 (Ajiki et al. 2002; Shioya et al. 2002). 
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(Kaifu 1998) at Mauna Kea Observatories. The Suprime-Cam 
consists of ten 2k x 4k CCD chips and provides a very wide field 
of view; 34' x 27'. 

In this survey 7 , we used the narrow-passband filter, A/B816, 
centered at 8150 A with the passband of AA(FWHM) = 120 
A; the central wavelength corresponds to a redshift of 0.24 for 
Ha emission. We also used broad-passband filters, B, Rc, Ic, 
and z' ■ All the observations were done under photometric con- 
dition and the seeing size was between 0.7 arcsec and 1 .3 arcsec 
during the run. 

The individual CCD data were reduced and combined us- 
ing IRAF and the mosaic-CCD data reduction software devel- 
oped by Yagi et al. (2002). Detail of flux calibration will be 
presented in Fujita et al. (2003b). The combined images for 
the individual bands were aligned and smoothed with Gaussian 
kernels to match their seeing sizes. The region contaminated by 
fringing pattern is masked. The final images cover a contiguous 
706 arcmin 2 area with a PSF FWHM of 1 ."4. A total volume of 
3.9 x 10 3 Mpc 3 is probed in ourA/B816 image. The net exposure 
times of the final images are 28, 80, 56, 86, and 600 minutes for 
B, Rc, Ic, z 1 , and M?816, respectively. The limiting magnitudes 
are B = 26.6, R c = 26.2, I c = 25.9, z' = 25.3, and A/B816= 26.0 
for a 3a detection on a 2". 8 diameter aperture. 

Catalogues of the objects are made using SExtractor (Bertin 
& Arnouts 1996). The objects are detected using the double- 
image mode: the narrowband frame is used as a reference im- 
age for detection and then the flux is summed up in 14-pixel 
diameter apertures in all images. This aperture size is 2". 8, i.e., 
2 x FWHM of the stellar objects in the final image of /?c-band. 



3. RESULTS 
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3.1. Selection ofNB816-Excess Objects 

Since the central wavelength of the Ic is bluer than that of the 
NB8 16 filter, we calculated magnitude that we refer to as the "Iz 
continuum", using a linear combination (Iz = 0.76/c+0.24z') of 
the I c and z' flux densities; a 3 a photometric limit of Iz ~ 26.0 
in a 2". 8 diameter aperture. This enables us to more precisely 
sample the continuum at the same effective wavelength as that 
of the M?816 filter. Following the manner of the previous sur- 
veys using narrowband filter and taking the scatter in the Iz- 
M?816 color and our survey depth into account, candidate line- 
emitting objects are selected with the criterion of 7z-M?816 > 
max(0.1, 3cr error of Iz-NBUb); note that Iz-NBUb = 0.1 cor- 
responds to EW bs ~ 12 A. We compute the 3a of the color 
as 3(7^-^816 = -2.51og(l - v / (/3a™ 8 J 2 + (/3 CT;z )7/A«8i6). This 
criterion is shown by the solid and dashed lines in Figure 1. 
There are 1224 sources that satisfy the above criterion. They 
are brighter than the limiting magnitude at each band. 

3.2. Selection ofNB816-Excess Objects at z~ 0.24 

A narrowband survey of emission line galaxies can poten- 
tially detect galaxies with different emission lines at different 
redshifts. If the source redshift and the rest frame wavelength of 
the line act to place it inside the narrowband filter, the line will 
be detected if it is sufficiently strong. The emission lines we 
would expect to detect are Ha, H/3, [O III] AA4959,5007, and [O 
II]A3727 (Tresse et al. 1999; Kennicutt 1992b) as the narrow- 
band filter passband is too wide to separate [N ll]AA6548,6584 
from Ha. In Table 1 we show the different redshift coverage 
for each line. 

In order to distinguish Ha emitters at z ~ 0.24 from 
emission-line objects at other redshifts, we investigate their 
broad-band color properties comparing observed colors of our 
1224 emitters with model ones that are estimated by using 
the population synthesis model GISSEL96 (Bruzual & Chariot 
1993). In Figure 2, we show the B-Rq vs. Rq-Ic color-color 
diagram of the 1224 sources and the loci of model galaxies. 
Then we find that Ha emitters at z ~ 0.24 can be selected by 
adopting the following three criteria (see for detail, Fujita et al. 
2003b); (1) B-R c > 2(R C - / C ) + 0.2, (2) z'-NB9,l6 > 0, and 
(3) / C -/VB816 < 0.75(z'-A«816) + 0.35. These criteria give 
us a sample of 348 Ha emitting galaxy candidates. In particu- 
lar, the first color criterion is a quite nice discriminator between 
Ha emitters at z ~ 0.24 and other emitters at different redshifts 
and thus it is expected that there is few contamination in our 
sample. 

We adopt the same method as that used by Pascual et al. 
(2001) to calculate Ha equivalent width. The flux density in 
each filter can be expressed as the sum of the line flux and the 
continuum flux density (the line is covered by both filters): 

/iw> = /c+^,and/ ft =/ c + 0.76-^ (1) 

with fc the continuum flux density; Fl the line flux; ANB and 
Al the narrow and I c band filter effective widths and f NB and 
fi z the flux density in each filter. Then the line flux, continuum 
flux density, and equivalent width can be expressed as follows: 

F L = ANB fm ~ f,z - , (2) 

L 1-0.76(AATB/A/)' v ' 

f Iz -0.76f NB (ANB/Al) 

\-Q.16(ANB/Al) 1 

^ = ANb( fNB ~ f ': - 

fc \f Iz -Q.16f NB (ANB/Al) 



fc = 



EW, 



obs 



(3) 



(4) 



where Al = 1269 A for the I c band and ANB = 120 A for the 
NB816 band. The minimum line flux in our 348 sources is 
1.2 x 10" 17 ergs s -1 and the minimum continuum flux density 
reaches 3.0 x 10" 30 ergs s" 1 Hz" 1 (w 25.2 in AB magnitude). 



3.3. Ha Luminosity 

In order to obtain the Ha luminosity for each source, we 
correct for the presence of [N II]AA6548,6584 lines. Further, 
we also apply a mean internal extinction correction to each ob- 
ject. For these two corrections, we have adopted the flux ratio 
of /(Ha)//([N II]AA6548,6584) = 2.3 (obtained by Kennicutt 
1992a; Gallego et al. 1997; used by TM98; Yan et al. 1999; 
Iwamuro et al. 2000) and A Hct = 1. 

We also apply a statistical correction (28%; the average value 
of flux decrease due to the filter transmission) to measured flux 
because the filter transmission function is not square in shape. 
The Ha flux is given by: 



/cor(Ha) = /(Ha+ [Nll])x 
x 10°' 4Ah 



/(Ha) 



/(Ha)+/([NII]) 
x 1.28. 



(5) 



We calculate /(Ha+ [ N II]) from the total magnitude; note that 
the aperture size is determined from the /c-band image. Finally 
the Ha luminosity is given by L(Ha) = 4-irdl/ cor (Ha) using the 
redshift of the line at the center of the filter z = 0.242 and the 
luminosity distance, <f L . 



4. LUMINOSITY FUNCTION AND STAR FORMATION RATE 

4.1. Luminosity function of Ha emitters 

In order to investigate the star formation activity in galax- 
ies at z ~ 0.24, we construct the Ha luminosity func- 
tion (LF) for our Ha emitter sample by using the relation, 
$(logL(Ha))AlogL(Ha) = £V 1 /V\ where V' is the comoving 
volume, and the sum is over galaxies with Ha luminosity within 
the interval [logL(Ha) ± 0.5AlogL(Ha)]. We take account of 
the filter shape in the computation of the volume. The correc- 
tion can be as large as 13% for the faint galaxies [L(Ha) < 10 41 
ergs s -1 ] as compared to that for the brightest galaxies. Figure 3 
shows the result with the Schechter function (Schechter 1976) 
fit to the z < 0.3 Ha LF measured by TM98 (which is charac- 
terized by a = -1.35, <j>* = 10" 2 56 Mpc" 3 , and L* = 10 4L92 ergs 
s" 1 ; note that we convert the parameters to those of our adopted 
cosmology). 

Fitting Ha LF of L(Ha) > 10 40 ergs s" 1 with a Schechter 
function, we find that the best-fitting parameters are a = 
-1 .53 ± 0.15, log cj>* = -2.62 ± 0.34, and log/.* = 41 .95 ± 0.25. 
We find that the faint end slope of our LF is steeper than that 
obtained by TM98. This difference may be attributed to the dif- 
ferent source selection procedures between TM98 and ours; i.e., 
TM98 used /-selected CFRS galaxies lying at redshift below 
0.3 with the limiting magnitude of Iab = 22.5, while we used 
A^?816-selected galaxies with the much deeper limiting magni- 
tude of Ic — 25. In order to check this possibility, we have made 
LFs of I c < 24, I c < 23, and / c < 22 sources (Figure 4), using 
our /c-selected catalog. It is shown from this analysis that the 
LF for galaxies with a brighter /c-band limiting magnitude tend 
to have a shallower faint end slope although the brightest end 
of LFs is not affected by the /c-band limiting magnitude. The 
difference between LF of TM98 and our data is originated from 



Ha luminosity functions at z ~ 0.24 
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this point. In addition to this, our survey is concentrated in a 
small redshift range around z ~ 0.24 while the redshift range 
of TM98 is z < 0.3. This may also affect the LF shape. It 
should be noted that our result is consistent with that obtained 
by Pascual et al. (2001) although their data are available only 
for brighter objects with log L(Ha) > 42. 

4.2. Luminosity density and star formation rate density 
The Ha luminosity density is obtained by integrating the LF: 

/•OO 

£(Ha)= / $(L)M = r(«+2)^L„ (6) 
Jo 

We then find a total Ha luminosity per unit comoving vol- 
ume 10 39 - 65 -" " ergs s" 1 Mpc" 3 at z « 0.24 from our best fit 
LF. The errors quoted here are the standard deviations taking 
into account that the three Schechter parameters are correlated. 
Our estimate is quite similar to that obtained by Pascual et al. 
(2001), 10 39 - 79±0 09 ergs s" 1 Mpc" 3 . 

Not all the Ha luminosity is produced by star formation, be- 
cause active galactic nuclei (AGNs) can also contribute to the 
luminosity. The amount of this contribution depends on the se- 
lection criteria. For example, it is 8-17% in the CFRS low-z 



sample (Tresse et al. 1996), 8% in the local UCM (Gallego et 
al. 1995), and 17-28% in the 15R survey (Carter et al. 2001) for 
the number of galaxies. Following the manner of Pascual et al. 
(2001), we adopt that AGNs contribute to 15% of the luminos- 
ity density. We correct this AGNs contribution and then obtain 
the corrected luminosity density, 10 39 58 - 17 ergs s _1 Mpc -3 . 

The star formation rate can be estimated from the Ha lumi- 
nosity using SFR = 7.9 x lO- 42 L(Ha)M yr _1 , where L(Ha) is 
in units of ergs s _1 (Kennicutt 1998). Thus, the Ha luminos- 
ity density can be translated into the SFR density of psfr — 
0.036!g;j]^M o yr 1 Mpc" 3 (with AGN correction 0.03 1^;^ 
M C3 yr -1 Mpc -3 ). Figure 5 shows the evolution of the SFR den- 
sity from z = to z = 2.0. The SFR density measured here is 
higher than that of TM98. This is mainly due to the difference 
of faint end slope as shown before. Our result is consistent with 
the strong increase in SFR density from z = to z — 1 . 

We would like to thank the Subaru Telescope staff for their 
invaluable help. We thank the referee, Jesus Gallego, for his 
detailed comments that improved this article. We also thank T. 
Hayashino for his invaluable help. This work was financially 
supported in part by the Ministry of Education, Culture, Sports, 
Science, and Technology (Nos. 10044052 and 10304013). 
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Table 1 

Emission lines potentially detected inside the narrowband. 



Line 


Redshift range 






V x 10 4c 




Zl < Z < Z2 




(Mpc) 


(Mpc 3 ) 


Ha 


0.233 0.251 


0.242 


1220 


0.39 


[O iii]AA4959,5007 


0.616 0.640 


0.628 


3740 


2.39 


H/3 


0.664 0.689 


0.677 


4100 


2.62 


[O n]A3727 


1.17 1.20 


1.19 


8190 


7.30 



a Mean redshift. 
b Luminosity distance. 
c Comoving volume. 




NB816 

FIG. 1 . — Objects detected down to WB816 = 26 in the ,/VB816-selected catalog. The horizontal solid line corresponds to Iz-NB&16 = 0.1. Dashed line shows the 
distribution of 3a error. Because bright sources are saturated in the continuum frame, Iz-NB&16 is not zero at NB816< 19. 



Ha luminosity functions at z ~ 0.24 




-0.2 0.2 0.4 0.6 0.8 1 1.2 

Rc-Ic 

FIG. 2. — Plot of B-Rq vs. i?c — /c for the 1224 sources found with emitter selection criteria (gray circles). Colors of model galaxies (SB, Irr, Scd, Sbc, and E) 
from z — to z = 1.2 are shown with blue lines. Colors of z = 0.24, 0.64, 0.68, and 1.18 (for Ha, [O III], H/3, and [O II] emitters, respectively) are shown with red, 
green, purple, and yellow lines. We select the sources above the black line as Ha emitters. 



SFR M sun yr 1 




39.5 40 40.5 41 41.5 42 42.5 43 
log[L(Ha) erg s _1 ] 

FIG. 3. — Luminosity function at z ~ 0.24 (open circles). Solid line is the best fitted Schechter function. The data points with L(Ha)> 10 40 ergs s~' are used in 
the LF fitting. The TM98 Ha luminosity function at z < 0.3 is shown with dashed line. The data points of Pascual et al. (2001) are also shown by filled triangles. 
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FIG. 4. — Effects of the threshold of 1q magnitude in our Ha luminosity function. The solid and dashed lines have the same meanings as those in Figure 3. LFs of 
/c-selected Ha emitting galaxies are shown with triangles (Iq < 24), circles (7c < 23), and squares (7c < 22). 
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FIG. 5. — Star formation rate density at z ~ 0.24 derived from in this study (filled circle; ref. - 6) shown together with the previous investigations compiled by 
Trentham et al. (1999) and Tresse et al. (2002). Filled symbols show SFR derived from Ha LFs obtained by spectroscopy (triangles) and narrowband imaging 
(circles). Open circles are SFR density derived from UV LFs which includes no dust correction. References: 1 - Gallego et al. (1995), 2 - Gronwall (1999), 3 - 
Sullivan et al. (2000), 4 - TM98, 5 - Pascual et al. (2000), 7 - Tresse et al. (2002), 8 - Hopkins et al. (2000), 9 - Yan et al. (1999), 10 - Moorwood et al. (2000), 
1 1 - Treyer et al. (1998), 12 - Lilly et al. (1996), 13 - Connolly et al. (1997) 



